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Abstract: The study of exact quasi- normal modes [QNMs], and their associated 
quasi-normal frequencies [QNFs], has had a long and convoluted history -- replete 



with many rediscoveries of previously known results. In this article we shall collect 



and survey a number of known analytic results, and develop several new analytic 
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results — specifically we shall provide several new QNF results and estimates, in a 

form amenable for comparison with the extant literature. Apart from their intrinsic 



interest, these exact and approximate results serve as a backdrop and a consistency 
check on ongoing efforts to find general model-independent estimates for QNFs, and 
general model-independent bounds on transmission probabilities. Our calculations 
also provide yet another physics application of the Lambert W function. These ideas 
have relevance to fields as diverse as black hole physics, (where they are related 
to the damped oscillations of astrophysical black holes, to greybody factors for the 
Hawking radiation, and to more speculative state-counting models for the Bekenstein 
entropy), to quantum field theory (where they are related to Casimir energies in 
unbounded systems), through to condensed matter physics, (where one may literally 
be interested in an electron tunelling through a physical barrier). 
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1. Introduction 



The investigation of exact quasi- normal modes [QNMs], and their associated quasi- 
normal frequencies [QNFs] , has had a long and convoluted history, replete with many 
rediscoveries of previously known results [|IJ. Physically there are strong connections 
to black hole physics, at least some versions of quantum gravity, and quantum field 
theory. There are also strong connections between QNMs, QNFs, and transmission 
resonances — which is why analytic understanding of QNFs is intimately related to 
knowledge of exact solutions for transmission amplitudes and transmission probabili- 
ties. In this article (focussing mainly on QNFs) we shall collect and survey a number 
of known analytic results, develop several new analytic results, and develop sev- 
eral new perturbative estimates, in a form amenable to comparison with the extant 
literature. 

In particular we shall first discuss standard textbook fare such as the delta- 
function potential, double-delta-function potential, and asymmetric double-delta- 
function potential [g |, g|, |, [7] §, g [TD|, 0, p], |l|, [jj, using them to illustrate 
transmission resonances, damped modes, and QNFs — with already some significant 
new results even at this elementary level. We shall develop exact analytic results 
for the QNFs of the double-delta-function potential in terms of the Lambert W 
function, and develop some perturbative estimates for the QNFs of the asymmetric 
double-delta-function potential. Secondly, we shall turn to the step barrier and the 
symmetric and asymmetric rectangular potential barriers 0, |3], f|, ^, |5], [7|, P|, |9|, |TU, 



TT| , |12| , p~3| , [14fl , again extracting some new results on QNFs and damped modes. 

Finally we shall then turn to something more challenging, the Eckart poten- 
tial |1| and its common simplifications — the tanh(-) and sech (•) potentials. The 
history of the Eckart potential is particularly complicated: Apart from the special 
case of the Morse potential (which pre-dates Eckart's work by one year Jl6[), other 
special cases and equivalent reformulations of Eckart's results (such as the Rosen- 
Poschl-Teller 
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Manning-Rosen | 19fl , Manning |[20| , Hulthen [21 1 , Ti 



etz |22|| , and Hua |23| potentials) post-date Eckart's work by anything ranging from 
several years to six decades. While the QNFs for the sech 2 (-) potential are quite 
standard and well-known, the QNFs for the tanh(-) potential and general Eckart 
potential are at best only implicitly given in the extant literature. 

Apart from their intrinsic interest, these exact results serve as a backdrop and 
a consistency check on ongoing efforts to find general model-independent bounds on 
transmission probabilities J25|, 26[ |27|, |28|, |29], |30fl , and general model independent 
estimates for QNFs [[H], [32], [33], [34], [35], [36], [37], |3~8fl . Thus these ideas have relevance 
to fields as diverse as black hole physics (where they are related to the QNFs con- 
trolling the damped oscillations of astrophysical black holes, to greybody factors 
for the Hawking radiation, and more speculatively to state-counting models for the 
Bekenstein entropy) , quantum field theory (where analysis of the QNFs is an impor- 
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tant technical step in calculating Casimir energies for unbounded systems |39], |40fl), 
through to condensed matter physics (where one may literally be interested in an 
electron tunelling through a physical barrier). 

We will be particularly interested in the QNFs (complex energies corresponding 
to purely outgoing waves in both spatial directions) [[35], [36], |37|, [3"8|, p5|, f£0|, |4T| , ^ , (4~3| , 
|4"4] , pE5|, P| [4"T| , [4"8| , and the closely related transmission resonances (real energies where 
the transmission probability is unity). Formally, QNFs are most easily found by 
looking for complex frequencies where the transmission amplitude becomes infinite. 

Experience obtained many independent fields (including black hole physics) has 
shown that it is quite common for the QNFs to be approximately of the form 



UJ r , 



(offset) + m(gap) + 0(1/ n) 



XI) 



where the "offset" is generally a complex number and the "gap" is typically a real 
number g|, S IS © S Ml U 0, I 



43, 44, 45, M 



We shall 

be particularly interested in checking for such asymptotic behaviour in the specific 
models we consider below. In particular the asymptotic behaviour of the QNFs for 
the double-delta potential is more subtle than one might have expected. 

2. Conventions 



To set the stage, we are interested in solving the time-independent Schrodinger 
equation 

/)-' tr- 
ip (x) = Eip(x). 



IiL^L + v(x ) 

2m dx 2 



In regions where the potential is zero the wavefunction takes the form 



ip(x) 



exp(±ikx) 
7^ 



2L.2 



E 



2m 



(2.1) 



(2.2) 



We shall define outgoing modes by 

exp(— ik\x\) 



ip(x) — > 



Vk 



oo. 



(2.3) 



QNFs are in this non-relativistic context more properly called quasi-normal wavenum- 
bers, and correspond to Im(k) > so that 



exp(Jm(/c)|x|) 
(x)\ -¥ — f1 ^ — ' -¥ oo; \x\ ->■ oo. 



(2.4) 



(Other sign and phase conventions are also in common use, and there is no universal 
agreement as to the "best" set of conventions, but all authors agree that the quasi- 
normal modes are growing at spatial infinity and so are non-normalizable.) Once 



-3- 



one has found the quasi-normal wavenumbers, the associated quasi-normal energies 
are (in this non-relativistic context) simply given by 



E, 



QNF 



h 2 k 2 QNF /(2m) 



(2.5) 



For symmetric situations where the potential has the same limit V(x) — > V^ at 
x — > ±00 the scattering problem is characterized by the asymptotic behaviour 

exp(— ik^x) 



t 



A (i 



[x — > +00] 



■ip(x) ->■ < 



exp(—ikoox) expiikoox) 

— \- r - 



(2.6) 



[x ->■ 



-00 . 



■y/k -y/k 

Should the potential have distinct limits V(x) — > ^±00 as x — > ±00 then one needs 
to distinguish the asymptotic wavenumbers k± 0O and the the scattering problem is 
characterized by the more complicated asymptotic behaviour 



exp(— ik A 



-,x) 



ip(x) ->■ < 



exp(—ik- 



x 



+00 



-,x) 



exp(ik. 



-,x\ 



\fk- 



\fki 



\x — > — 00 . 



(2.7) 



If one wishes to work directly in a relativistic context then instead of the Schrodinger 
equation the relevant PDE is the closely related 



+ 



+ V(x 



dt 2 dx 2 
For relativistic wavefunctions of the form ifj(t,x 

d 2 



dx' 



+ V(x) 



ij>(t,x) = 0. (2.8) 

1 = exp(iujt)i/)(x) this reduces to 
ip{x) = u 2 ip(x). (2.9) 



This is formally equivalent to the situation for the non-relativistic time-independent 
Schrodinger equation under the substitutions h 2 /(2m) — > 1 and E —$■ u 2 . Thus there 
is no particular need to treat the relativistic situation separately. We shall phrase 
the discussion below in terms of the non-relativistic problem with the understand- 
ing that the fully relativistic situation can easily be recovered with the appropriate 
substitutions. 

3. Selected potentials leading to exact results 

We start by looking at some extensions of textbook results and connecting them 
back to the commonly conjectured "(offset) -Mn(gap)" behaviour for highly damped 
QNFs. 
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3.1 Delta- function potential 

For a delta function potential take 

V(x) = a5(x). (3.1) 

Since the potential is zero for x 7^ we can in this region relate the energy to the 
wavenumber via 

h 2 k 2 

E = . (3.2) 

2m v ; 

It is extremely useful to define 

h = -p-. (3.3) 

In this case the transmission amplitude is well known to be (see for instance Baym || 
or Gasiorowicz ||) 

k 
where this expression holds even for complex k. For real k (and hence real E) the 
transmission probability is (see for instance 0, [|) 

T=^ F . (3.5) 

k 2 

There are no transmission resonances for this potential. (Speaking rather loosely one 
could formally view the limit k — y 00 as a transmission resonance, as T — y 1 in this 
limit — but this is not standard nomenclature.) 

The QNFs are located by finding the poles in the transmission amplitude t. 
There is only one QNF, which is pure imaginary, (and so strictly speaking should be 
called a damped mode rather than a QNF). It is given by 

h 2 k 2 

&qnf = iko) E QNF = — —A (3.6) 

2m 

For a > the delta function is repulsive, and Im(k) > for the formal QNF, 
corresponding to a damped-mode QNF. On the other hand, for a < the delta 
function is attractive — then Im(k) < for the formal QNF, so the formal QNF is 
actually a bound state. Because the width of the delta function is zero, the "gap" 
is infinite, and the "higher" QNFs are effectively driven off to imaginary infinity - 
only the "lowest" QNF survives. 

3.2 Double-delta-function potential 

For the double delta function 

V(x) = a{5(x-a) +5(x + a)}, (3.7) 
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the transmission amplitude is known to be (see for instance Galindo and Pascual 

A- 2 



t 



(k — iko) 2 + &Q exp(— 4ika) 
even for complex k. Here we again set 



At: 



ma 



For real k the transmission probability is (see for instance 

T- l 



Equivalent ly 



T 




(3. 



(3.9) 



(3.10) 



k 
1 + 4 -| [k cos{2ka) + k sin(2fca)] 2 



(3.11) 



The transmission resonances are located by 

T = 1 <^=^> k = —k ta.n(2ka) 



(3.12) 



This is a transcendental equation with an infinite family of exact solutions. For large 
real k approximate solutions are 



The QNFs are located by 

t = oo 
so that 

Alternatively 



2ka ps ( n + - 1 7T. 



(A; — zA;o) + k exp(-Aika) = 0, 



(3.13) 



exp(4iA;a) 



ik 
k Q 



exp(-2ika) =■■ ±( 1 + — ). 



/A' 



(3.14) 



(3.15) 



(3.16) 



One obvious formal solution is k = 0, but this is physically uninteresting and does 
not correspond to a true physical QNF. The general solution to the QNF condition 

is 

W{±2k a e 2koa ) 



k 



QNF 



% { k 



2a 



(3.17) 
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Here W(x) is the Lambert W function implicitly defined by W[x) e w ^ = x. (For 
a general discussion of the Lambert W function see [|9|, [5^, [51], ^2|, |53], ^4|.) Now 
the Lambert W function is, like the logarithm or many other complex functions, a 
multi- valued function depending on which particular branch one is dealing with. The 
real branches are Wo{x) for x G (— e _1 , oo) and W-i(x) for x G (— e _1 , 0). For these 
two branches W{o,_i} (xe x ) = x, while on the other hand W^{o,-i}( — xe - x ) is nontrivial. 
This now leads to one trivial QNF at k = 0, plus a second nontrivial QNF at 

&QNF = « S Ko g - > • (3.18) 

This particular QNF is pure imaginary (either a damped mode or a bound state) 
whenever the W function is real, that is whenever 

2ak < W(e- 1 ) = 0.2785... (3.19) 

This particular exact result for /cqnf is already a rather nontrivial analytic result 
that does not appear to be well-known. As a — > 0, so that the 2 delta functions 
merge and W(x) = x + 0(x 2 ), this reproduces the previous single delta function 
result. In addition, from the other (guaranteed complex) branches of the W function 
one obtains an infinite tower of QNFs — specifically 

KqnfW = %< « > . (3.20) 

The corresponding quasi-normal energies are 

p n ^ It iy n (±2fc qe 2 ^) ) 2 

^qnfW = -— |fco ^ ) • (3.21) 

Thus already in this simple situation, as soon as one has non-zero separation a 
between the two delta-functions, then an infinite tower of QNFs arises. As the 
"width" a shrinks to zero all but the lowest lying QNF are driven off to infinity. 
Quantitative asymptotic information can be extracted by first noting 

W n (x) = -lnW n {x) + lnx + i27m. (3.22) 

By recursively iterating this expression one can derive the Comtet expansion (see for 
example [5(|). The leading terms of this Comtet expansion can best be recast as 



WJx) = Inx + %2-nn - ln(lnx + %2-nn) + O J , x + z — I , (3.23) 

^ mx + i2im J 

where the remainder term slowly goes to zero as \n\ —¥ oo. Let s = ^±|g{0,1} 
depending on the sign ±. Then 

W n (±2k a e 2koa ) = 2k a + \n(2k a) + i(2n + s)vr 

+ In [2k a + ln(2k a) + i(2n + s)n] + ... (3.24) 
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with slowly decreasing remainder as \n\ — > oo. Consequently 
A; QNF ( n ) = (2n + s)7r _. ln(2fc a) _ % [n ^ + ln(2A;oa) + ?(2n + s)?r] + _ _ _ (3-25) 

Note that it is in this case the real part of the QNF that is equi-spaced as as \n\ — \ oo, 
not the imaginary part. Note further that the sub-leading term is logarithmic, not 
0(l/n). This indicates that our intuition regarding the (offset) + m(gap) behaviour 
built up from black hole mechanics is not quite as universal as one might hope — at 
the very least we will have to permit imaginary gaps and offsets, and more general 
sub-leading terms. 

3.3 Asymmetric double-delta-function potential 

For the asymmetric double delta function 

V(x) = a_ 5(x — a) + a + S(x + a), (3.26) 

with the definitions 

ma+ , ma_ 

k + = -^±; fc_ = —■ (3.27) 

we have the transmission amplitude 

P 
^ _ /o 2 g\ 

(k-ik + )(k-ik_) + k + k-e- 4ika ' 
As a — > one has 

' "» S-+IT- < 3 ' 29 > 

in agreement with the single delta function case. Also, for k + = k_ — > fc this result 
reproduces the amplitude for the symmetric double-delta potential. For real k the 
transmission probability T can be written as 

tf 

k 4 + k 2 {k\ + k 2 _ ) + 2k%kl + 2k + k_[(k 2 - k + k_) cos(4ka) + k(k + + fc_) sin(4A;a)] ' 

(3.30) 

This can be recast as 

f 

T = 

ffc_L — kJ) 2 Ak±k_ 

1 + v — — - H i— [k cos(2ka) + k + sin(2£;a)] [k cos(2ka) + £;_ sin(2£;a)] 

k 2 k 4 

(3.31) 
As a — > one has 

T ^ k* + Z+k_ r (3 ' 32) 



in agreement with the single delta function case. Also, for k + = £;_ — > ko this result 
reproduces the amplitude for the symmetric double-delta potential. 

For k + 7^ k_ (and real energies) there are no true transmission resonances as T 
never quite reaches unity. However approximate transmission resonances, where one 
has a local maximum T < 1, are located at the local minima of the quantity 



(fc+ - fc- 
Aklkl 



k 

cos(2ka) H sin(2A;a) 

k 



k_ 

cos(2ka) H sin(2ka) 

k 



> 



0. (3.33) 



For large k the location of these approximate transmission resonances approximates 
to 

cos(2A;a) « 0. (3.34) 

That is, there is a family of approximate transmission resonances given by 

2ka^ [n+-j n. (3.35) 

Near these approximate transmission resonances 

T W (k + - k.f Aklk 2 ■ (3 ' 36) 



A; 2 A; 4 

So these approximate transmission resonances become asymptotically exact as the 
momentum becomes large. 

Turning to complex energies, the QNFs are located by 

(k - ik+)(k -ik-) + k + k^ e - 4ika = 0. (3.37) 

One formal solution is k — 0, but this does not correspond to a true physical QNF. 
The non-trivial QNFs do not seem to be explicitly calculable even in terms of the 
Lambert W function, but we can develop a perturbative result for k + ~ k_. 

To start the discussion we can at least see that one of the non-trivial QNFs is 
pure imaginary, (that is, a bound state/purely damped mode). Writing k — i\k\ we 
have 

(|ib| - fc+)(|A;| - fc_) = k + k_e~ ma . (3.38) 

Graphically, it is easy to see that for k + ^ k_ there will always be one trivial solution 
at k — 0, plus one distinct non-trivial QNF with \k\ > max{/c + , A;_}. 
More generally, the exact QNF condition can be rewritten as 

{2ia{k - ik + )e 2ika ) x {2ia{k - ik^)e 2ika } = 4a 2 k + k., (3.39) 



or 



{2ia{k - ik + )e 2i{k ~ ik+)a ) x {2ia{k - ik„)e 2i{k - ik - )a ) = Aa 2 k + k_e 2{k++k - )a . (3.40) 
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Defining 



C = 2a^kle {k++k - )a , (3.41) 

we have (again in terms of the Lambert W function) 

k = ik + -i w ^ c r +e) = ik.-i w ^ e - e \ (3.42) 

where we have the consistency condition that 

2a(k + - k-) = W n (C e +e ) - W n (C e" 6 ), (3.43) 

and for emphasis we have used W n to indicate that there is an infinite collection of 
QNFs hiding in the various branches of the Lambert W function. For k + ~ A;_ one 
has a perturbative result 

a(k + -k^ = W' a (C )Q + O(Q 3 ), (3.44) 

whence 



This easily leads to the self-consistent estimate 

, \ k + + k_ W n (C ) 2 \ 

With a bit more work one can extract the next higher-order estimate 



(3.46) 



, / x .\k + + k^ W n (C ) a{k + -k_) 2 4 . 

W) = .|— jj- " 4 , y „ (Co)|1 + w „ (Co) ] + 0(1*+ - *-l > I 

(3.47) 

These results for the perturbative estimation of an infinite tower of QNFs appear to 
be new. The corresponding quasi-normal energies are 

P / n & f k + + k - W n {C ) a{k + - A;_) 2 4 ^ 

^ QNF(n) = -^{—2 2^~ ~ AW n (C )[l + W n (C )] + °^ - k ~ ] ) 

(3.48) 

We could again use the dominant terms of the Comtet expansion, now for W n (Co), 
to determine the asymptotic behaviour as \n\ — > oo. We suppress the details as 
the overall behaviour is qualitatively similar to that for the symmetric double-delta 
potential considered above. 

In contrast, for small a one has a different perturbative relation (now for the 
lowest QNF) 

A; QNF = i(k+ + k-) - 4ik + k_a + 0{a 2 ), (3.49) 
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with 



,2 



E QNF = -^-{(k + + k^)~ 4k + k_a + 0(a 2 )} 2 , (3.50) 

which is compatible with the result for a single delta-function. 

In summary, all the explicit models we have seen based on one or two delta func- 
tion potentials lead to a solitary isolated damped mode/bound state (pure imaginary 
/cqnf), and in addition the situation with two separated delta function potentials 
leads to an infinite tower of generally complex QNFs. 

3.4 Single-step potential 

The single step potential has the form (see for example Landau-Lifshitz ||): 

(for x > 0): 
V(x) = { (3.51) 

(for x < 0). 

Let us write 

e= 2j^E. ,. J^-".) . 3 = ^ = *"-«"■ (3.52) 

Then (see for example Landau-Lifshitz ||, Messiah 0, or Merzbacher ||, though 
note different flux conventions) 



* = g*. (3.53) 

k + q V ' 

If k and q are both real (so that E and E — Vq are both real and positive, which 
is the minimum requirement for a true scattering situation) then the transmission 
probability is 

T = w+w = l ~ w+w (3 ' 54) 

There are no transmission resonances for this potential. (Except formally in the limit 
E —¥ oo.) There are no QNFs for this potential. From the point of view of QNFs a 
single-step potential is uninteresting - - but this will change radically once two-step 
potentials are considered. 

3.5 Rectangular barrier 

The rectangular barrier (sometimes called the square well) has the form (see for 
example Landau-Lifshitz || or Schiff 0: 

(for \x\ < a): 
V{x) = { (3.55) 

(otherwise). 
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Where in general Vq can be either positive or negative. Let us write 

2 2mE 2 2m(E -V Q ) 2 2mV o ,22 / Q k R \ 

fc = -^-; g = ^ ; k = —^- = k -q . (3.56) 

The exact transmission amplitude is (see for example Brandsen and Joachain [ |T(i|j , 
or Messiah [0], though note change in phase conventions) 

Akq exp(2ika) 



{k + q) 2 exp(2zga) — {k — q) 2 exp(— 2iqa) 
For real k (real E) the exact transmission coefficient is 



(3.57) 



k 2 q 2 
T = k 2 q 2 + \(k 2 -q 2 ) 2 sm 2 (2qay (3 ' 58) 

which can be rewritten as 

T= g ~ Vo) (359) 

E(E - Vq) + \V 2 sm 2 (2^2m(E - V )a/h) 

For more details see (for example) the texts by Landau and Lifshitz [BJ , or Schiff M . 
We can re-write this as 

T = ; . (3.60) 

mVja 2 sm 2 {2^/2m(E-V )a/h) 

1+ 2EK 2 2m(E - V )a 2 /h 2 

The transmission resonances are defined by 

r = i ^ ,= -. (3.61) 

This now is the first time we see "families" of exactly solvable transmission resonances 
arising. In terms of the incident energy the transmission resonances occur at 

h 2 n 2 7r 2 

E = V + -S— 2"- 3 - 62 

sma z 
In contrast the QNFs are located by looking at complex wavenumbers and are defined 
by t = 00 corresponding to 

(k + q) 2 exp(2iqa) — (k — q) 2 exp(—2iqa), (3.63) 

that is 

(k + q) exp(iga) = ±(k — q) exp(-iqa), (3.64) 

which can be simplified to 

k = —iq tan(ga), or k = iq cot(ga), (3.65) 

and thence to 

ko = — iq sec(ga), or k = iq csc(qa) , (3.66) 

whence 

q = ik cos(ga), or q = —iko sin(ga). (3.67) 

There are two cases of interest for dealing with these transcendental equations. 
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3.5.1 Attractive potential 

For an attractive potential Vq < 0, and then /c = i\ko\ is imaginary. So for an 
attractive potential we would be interested in 

q = —\ko\ cos(ga), or q — \ko\ sin(ga). (3.68) 

Apart form the trivial solution at q — 0, these equations (sometimes) have solutions 
for real values of q, but always with \q\ < \k \. Since k 2 = k^+q 2 = — \ko\ 2 +q 2 < this 
implies pure imaginary values of k, which correspond to bound states (not QNFs). 

In addition, suppose Im{q) ^> 0, then cos(ga) ~ exp(— iqa)/2 and sin(ga) p» 
— exp(— iqa)/2. In either case the QNF condition approximates to 

\k I 
q ~ — exp(— iqa) (3.69) 



with approximate solution 

. W n (-i\ko\a/2) 
9qnf(w) « -i . (3.70) 

Hi 

This approximation should become increasingly accurate for higher values of Im(q), 
that is, for higher values of n. 

3.5.2 Repulsive potential 

On the other hand, for a repulsive potential Vq > 0, and then fc > is real. Let us 
first consider the situation q = i\q\ (that is, a purely damped mode/ bound state). 
Then we must solve 

\q\ = ko cosh(|g|a), or q = 0. (3-71) 

The zero mode is not a true QNF, and the physical QNFs are defined by 

\q\a = k a cosh(|g|a). (3.72) 

Depending on the precise value of k^a there will be two, one, or zero QNFs. For 
small koa there are two QNFs and one can perturbatively estimate the lower of these 
two QNFs by 

f 1 13 1 

g Q NF = ik 1 1+ ^ (M 2 + ^(M 4 + O([k af) j . (3.73) 

This corresponds to 



%\k = /7 u ,(/mk0 \ 1 + g(M 2 + g(M 4 + O{\k af) \ . 



(3.74) 
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The quasi-normal energy is 
-Eqnf = 



ft 2 k 2 f 4 92 

^(M 2 1 1 + |(M 2 + ||(M 4 + <?([M ( 



(3.75) 



In contrast, the higher of these two QNFs seems to have no simple perturbative 
expansion or other representation. For k a ~ 0.663 the two QNFs merge, (at qa ~ 
1.2), and for k a > 0.663 there are no longer any QNFs. 

On the other hand, for Im(q) ^> 0, we can again approximate the trigonometric 
functions by exponentials. The QNF condition approximates to 



iko , . , 
q m — — exp(— iqa) 



with approximate solution 



?QNf(™J 



.W n (k a/2) 



(3.76) 



(3.77) 



This approximation should become increasingly accurate for higher values of Im(q), 
that is, for higher values of n. 

3.6 Asymmetric rectangular barrier 

For the asymmetric rectangular barrier (sometimes called the asymmetric square 
barrier) 

{Vi, x < —a; 

V 2 , xe(-a,+a); (3.78) 

V3, x > a. 

We now define 

kt = y/2m(E-Vi)/h, 16(1,2, 3}, (3.79) 

where the ki are either pure real or pure imaginary, and construct a wave-function of 
the form (see for instance Messiah, though note change in normalization and phase 
conventions M) 



{exp(— ik\x) +rexp(ikix)} 



ip(x) 



{p exp(—ik2x) + q exp(ik2x)} 
texp(— ik^x) 

7e 



x < —a; 



x G (—a, +a): 



x > a. 



(3.80) 



The continuity conditions at points —a and +a give the values of r, p, q, and t. 
Without entering into the specific details of the calculation, we simply give the 
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probability amplitude t. A brief computation yields: 



2k2Vhki exp(i[fci + k 3 ]a) 

^2(^3 + fci)cos(2A; 2 a) + i{k\ + feiA;3)sin(2A;2a) ' 



This is more usefully represented as 



4£; 2 vM^3 exp(i[kt + fc 3 ]a) 



(A* + fc 2 )(fc 3 + /c 2 )e 2ifc 2« - (fci - fc 2 )(fc 3 - k 2 )e~ 2ik ^' 



(3.82) 



T sin(-)^0 = 71 , ,.% 2 = ^stepCfcl, fc 3 ). (3-84) 



When fci and fc 3 are both real, which is the minimum requirement to have a true 
scattering problem, the transmission probability is (see for example 0): 

4fci/c 2 fc 3 
= (£4 + k 3 ) 2 k 2 2 + [k 2 k 2 3 + k\{k\ -k\- k 2 )] sin 2 (2fc 2 a) ' (3 ' 83) 

The definition of a transmission resonance is now more subtle -- clearly something 
special happens when the sin(-) — y 0. Then T — y T step — at this point the transmis- 
sion probability for the asymmetric rectangular well reduces to that of the transmis- 
sion probability for a step potential with the same asymptotic behaviour at spatial 

infinity. 

4fcifc 3 

(h + h 

(In fact it is known that this step barrier transmission coefficient is a rigorous upper 
bound on the general transmission probability |25| .) We now define the closest we 
can get to a transmission resonance by asking 

T/T step = 1 -<=>- 2A; 2 a = nn. (3.85) 

In terms of the energy these "close as possible to transmission" resonances occur at 

ft n ir 
E = V 2 + t ^ T . (3.86) 

The QNFs are now determined by locating complex wavenumbers for which t — 00, 
corresponding to 

^2(^3 + k 1 )cos(2k 2 a) = — i(k 2 + kik 3 )sm{2k 2 a) . (3.87) 

That is 

/ , k 2 (k-i + k\) , 

tan(2 M = + i ^±^. (3.88) 

Let us now define constants k\ 2 and /c 23 (which are either pure real or pure imaginary) 
by 



ii2 = v 'M^) ; ^M^ (3 . 89) 
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so that 

h = sjh* + k\ 2 - h = y^i + k 2 23 . (3.90) 

Then the QNF condition becomes 

k 2 (y/k% + k\ 2 + y/kl + k 2 32 ) 

tan{2k 2 a) = +i -r^ \ ■ (3-91) 

( k 2 + v k 2 + k 12 y k 2 -\- k 32 j 

This implicitly determines k 2 (a, ki 2 , k 23 ) as a function of the parameters a, ki 2 , and 
k 23 . This can now be solved perturbatively for k 2 as a function of a, though the 
analysis is now considerably more delicate because there are several energy scales in 
play. Provided that both \k 2 \a <C 1, (so that we can safely approximate the tangent 
function by a straight line), and \k\ 2 — k 23 \/\k 2 2 + k 23 \ <C 1, (so that the potential is 
not too asymmetric), one can derive the estimate 

i (U2 _ p \2 9 p p 

(t, \2 _ l VH2 ^23) ^12^23 7,2 7,2 n 2 , n( 4\ (Q QO\ 

{k 2 ,QNF) ~ -7-3 7p ,p» " p r p ~ fc 12 fc 23 Q + ^l Q J- (3.92) 

^ u 1^12 "r ^23/ ^12 1" "-23 

Note that for small a this guarantees a pure imaginary QNF (damped mode/bound 
state). Insofar as the calculations overlap, this agrees with the perturbative QNF 
estimate for the symmetric rectangular barrier. The corresponding quasi-normal 
energy is 

h 2 (V.-V,) 2 2(V 2 -V 1 )(V 2 -V 3 ) 



EqNF V * 8ma 2 (2V 2 -V 1 -V 3 ) 2 (2V 2 -V 1 -V 3 

{v 2 -v x )(y 2 -v z ) , „. 4 



*/(W) + ° {a) - (3 ' 93) 

If we now want to move beyond this (small a) perturbative approximation, the alge- 
bra unfortunately appears to become intractable. We have not been able to extract 
anything useful in this case. 

3.7 Tanh potential 

Consider a smoothed step function of the form 

V {x) = V -°° I V+ °° + V+Q ° ~ V -°° tanh (-\ . (3.94) 



Define 



y 2m{h v±oo) r _ k+00 + k-oo ^ k +00 — re_oo , , 

k±oc — r ; k — - ; a — - . (^o.yoj 



The transmission amplitude is known analytically to be (see, for instance, || [14 

k T(ika) 2 

y/k +00 fc.oo r(zA; +00 a)r(iA;_ 00 o) ' 
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Furthermore, the transmission probability is known analytically to be (see for in- 
stance |, 0): 



sinh(7rA; +0O a) sinh(7r/c_ 00 a) 
sinh 2 (7rA;a) 



sinh (7rA a) 
sinh 2 (7r/ca) 



(3.97) 



There are no transmission resonances for this potential. (Except formally in the 
limit k — > oo.) To find the QNFs note that t — oo when the Gamma function in the 
numerator has a pole, that is when 



ika = —n, n G {0, 1, 2, 3, . . . }, 



(3.98) 



that is, when 

ia^-oo + k +oa )/2 = -n. (3.99) 

But we can eliminate one of the asymptotic wavenumbers in terms of the other plus 
the asymptotic potential difference V+oo — ^_oo. For instance 



k+oo + \Jk 2 +OQ + 2m(V +oc - V-oo)/h 2 = 2in/a , 
so that in terms of the wavenumbers at x — > +oo the exact QNFs are 



k 



QNF,+oo 



n) 



m(V +OD - l / _ 00 )g n 
2h 2 n a 



n > 



(3.100) 



(3.101) 



We can completely equivalently specify these exact QNFs in terms of the wavenum- 
bers at x — > — oo, in which case 



k 



QNF,-oo 



n 



m{y_ 



oo — > / +oo) a ' 



2h 2 n 



n 

+ - 

a 



n > 0. 



(3.102) 



Though the step from transmission amplitude to QNF is in principle straightforward, 
an explicit statement as to the exact location of these QNFs seems impossible to find 
in the extant literature. 

Note the asymptotic spacing as n — > oo: 



Ti 

&QNF,±ooW ^ i -■ 

a 



(3.103) 



Also note that as a — > all the QNFs are driven to imaginary infinity - - this is 
compatible with the behaviour of the single-step potential for which there are no 
QNFs. If one prefers to work with the non-relativistic quasi-normal energies then 



E, 



h 2 



V +00 + V_ c 



m 2 a 2 (V^ 



, " 2 , ^+oo ~r v -co '"■ " V+oo v-oo) J- / Q 1n/j \ 

QNF - = - W n + 2 8^ n 2 ' (3 - 104) 

Again, though the step from transmission amplitude to this (non-relativistic) quasi- 
normal energy is in principle straightforward, an explicit statement as to the exact 
location of these quasi-normal energies seems impossible to find in the extant litera- 
ture. 
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3.8 Sech 2 potential 

Consider a sech 2 potential of the form 

V(x) = V sech 2 (x/a) . (3.105) 

The transmission amplitude is known analytically to be (see for example Landau- 
Lifshitz § or Beyer ||): 



T (ika + \ + J\ - 2mV a 2 /h 2 ) T (ika + \ - J\ - 2mV a 2 /h 2 
1 = T(l + ika) f(ik) "' (3,1 ° 6) 

Furthermore, the transmission probability is known analytically to be (see for exam- 
ple Landau-Lifshitz fpfl): 



T = smh 2 [irV2mEa/h} 

sinh 2 [7r V2mEa/h] + cos 2 [±w^l - 8mV a 2 /h 2 } ' 



Alternatively 

sinh [nka] 



T 



sinh 2 [7rfca] + cos 2 ffyj — 2mVoa 2 /h 2 



(3.108) 



Transmission resonances occur but they are not now functions of energy or momen- 
tum. Instead T — > 1 for 



x/1 - 8mV a 2 /H 2 = 2n + 1. (3.109) 

That is, the potential is reflectionless when the depth of the well takes on one of the 
critical values 

The QNFs are in this case quite standard and are presented for instance in Beyer fl8| 
and many other places in the literature. We have 



\± J\-2mV a 2 /h 2 + n 
fc QN F,„ = *<( ) ne {0,1,2,3,...}. (3.111) 

Note the (offset) + m(gap) behaviour. Note that as a — > for Vo fixed all but one 
of the QNFs are driven to imaginary infinity (the one remaining QNF is driven to 
the unphysical value of zero). In contrast as a — > for VqCl held fixed all but one 
of the QNFs are driven to imaginary infinity and the one remaining QNF is driven 
to the unique QNF for a delta-function potential. The corresponding quasi-normal 
energies are 

h 2 [\ IT 



2ma 2 2 V 4 



£ Q NF,n = -t: ^ 1 o ± V 7 - 2mV a 2 /h 2 + n } . (3.112) 
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3.9 Eckart/ Rosen— Morse/ Morse— Feshbach potential 

We shall soon see that the Eckart potential can best be viewed as a linear combination 
of the (tanh) and (sech) 2 potentials. To set the stage we emphasize that many of the 
apparently different potentials commonly encountered in the literature are actually 
the same quantity in disguise. To start with, consider the following three potentials: 



V(x) = - T -^ 7 - n \ 2 ; £=-exp(2x/a). (3.113) 



Eckart (1930): 

i-e (i-0 5 

Rosen-Morse (1932): 

V(x) = A + Btanh(x/d) + C sech 2 (x / d) . (3.114) 

Morse-Feshbach (1953): 

V(x) = V cosh 2 /i {tanh([x — /jlL]/L) + tanh/i} 2 . (3.115) 

All three of these potentials are actually identical. To see this note that: 
4£ 4 4 1 



(1-0 2 "" it 1 ' 2 + i +1/2 ) 2 " [e-*/- + e*/°] 2 " cosh 2 (x/a) 
Similarly 



sech 2 (x/a). (3.116) 



2f l _j_ t 1 _ g2x/a g-x/a _ e x/a 

1 + 1^ = T3^ = 1 + e 2x/a = e -x/a + e x/a = -taBh(x/o). (3.117) 

This is enough to show 

(Eckart) <^=^ (Rosen-Morse). 



In fact in the Rosen-Morse article |17|], they cite Eckart |f[5[j , and describe Eckart's 



potential as begining "somewhat like" their own, but without noticing that the two 
potentials are in fact identical up to trivial redefinitions of the parameters. 

Now, for the Morse-Feshbach potential || , note that by a trivial shift of origin, 
x — > x + fiL, we have 

V(x) -> "i/ cosh 2 /i{tanh(x/L) + tanh/i} 2 , (3.118) 

which we can without loss of generality relabel as 

V{x) -»- V 1 {tanh(x/L) + D} 2 , 

= V x {tanh 2 (x/L) + 2Dtanh(x/L) + D 2 }, 

= Vi {-sech 2 (x/L) + 2Dt&nh(x/L) + D 2 + 1}, 

= V 2 sech 2 (x/L) + V 3 t&nh(x/L) + V 4 . (3.119) 
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This is enough to show 

(Morse-Feshbach) •<=>- (Rosen-Morse), 

and so all three potentials are completely identical up to trivial relabeling of the 
parameters and a shift in the zero of energy. 

In fact, including the offset, all three of these potentials can be written in any 
one of the four general forms below: 

V(x) = A + Btanh(x/a + 9) + Ctanh 2 (a;/a + 9), 
= A + [5c + C tanh(x/a + 9)f , 
B\ + C\ tanh(x/a) 



A + 
A + 



B 2 + Ci tanh(a:/a 

E x + Fi exp(-2x/a) 

_E 2 + F 2 exp(-2x/a) 



(3.120) 



Note that there is some redundancy here, but it is a useful redundancy. It makes 
it clear that the Eckart/ Rosen-Morse/ Morse-Feshbach potential is generally the 
square of a Mobius function, either of the variable tanh(x/o) or of the variable 
exp(— Ixl a). Thus implies that without loss of generality we can set either B\C 2 — 
C\B 2 or E\F 2 — E 2 F\ to some convenient constant (often unity). As long as E 2 and 
F 2 do not have opposite signs then the (Mobius) 2 potential does not exhibit any 
poles, and as long as all of Ei, -F\, E 2 , and F 2 are in addition nonzero it has the 
appropriate asymptotic behaviour to define a scattering problem. Other cases are 
more interesting as model potentials defined on a proper subset of the real line, and 
possibly with bound states. Sometimes one is interested in formally replacing a — >■ ia 
(and taking appropriate real parts) to obtain potentials based in combinations of 
ordinary trigonometric functions. Sometimes one is interested in formally replacing 
a — > —a, which can also lead to potentials with poles that are defined only on subsets 
of the real line. For our purposes we will only be interested in those specific forms 
of the (Mobius) 2 potential that lead to a well-defined scattering problem. 

To explicitly display the transmission amplitude and probability, and the trans- 
mission resonances and QNFs, it is convenient to settle on the standard notation 

V{x) = + Un^-j + —^- J - y (3.121) 

That is, the general Eckart potential is simply a linear combination "(constant) + 
(tanh) + (sech) 2 ". Now define the quantities 



k±00 = vgjjjgEEB; k = fc+oo+*-~ (3 . 122) 



h 
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The transmission amplitude is then known to be (see for example Eckart [15| or 
Morse-Feshbach [|J) 



V ^+00 fC- 



x 



r (ika + \ + J\ - 2mV a 2 /hA T (ika + ± - J\ - 2mV a 2 /h 2 



T(ik +QO a 



r(ik-c^a) 



(3.123) 



The transmission coefficient is |§ 

sinh( irk^oo a) sinh(7rfc +00 a) 



T 



sinh 2 (7rfca) + cos 2 



71 



1 _ 2mV a 2 

4 h 2 



(3.124) 



Note that this has appropriate limits as Vq — > where it reproduces the tanh poten- 
tial, and as V^^ — > V +00 where it reproduces the sech 2 potential. The closest one 
now gets to a transmission resonance is that T — > T tan h whenever the cos[-] — > in 
the above, that is, whenever the coefficient Vq of the sech 2 part of the potential takes 
on the special values 

h 2 



V 



2ma 2 



n(n + 1 



(3.125) 



Regarding the QNFs, we see that t — > 00 when 



ika H — ± \ 2mVoa 2 /h 2 = —n; 



ne {0,1,2,3,...}. (3.126) 



That is 



whence 



i(/c_oo + k +oc )a = ±\l 1 -^ (2n + 1) 



h 2 



\rv—OD ~T ^C^ 



1 / 8mV a 2 .2n + l 
a V h 2 a 



(3.127) 



(3.128) 



We can now rearrange this solely in terms of k +OCl by writing 



/, + x x /k 2 +00 + 2m(V +QO - V^)/h 2 = ±i-\l - ^f- + i^±l, (3.129) 

a V nr a 



implying 



k +00 + Jk 2 +oc + 2m(V +oc - V^/h 2 



'2n + 1 ± ^/l-8mV a 2 /W 



. (3.130) 
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Note that this has appropriate limits for the tanh and sech 2 potentials. Finally, 
because this is a simple quadratic equation, we can solve for the QNFs k(n) to 
obtain the exact QNFs 



, , , m(V +00 - V-^a/h 2 ( (2n + 1) ± ^1 - 8mV ayH 2 

^ONF +00 W = i ; H 

- ' l (2n+l)±^/l-8mV ayW 2a 

(3.131) 

A similar analysis in terms of the wavenumber k-oo at x — > — oo casts the exact 
QNFs in the form 



, m(^_oo - V +00 )a/h 2 (2n + 1) ± ^/l - 8mV a 2 /h 2 

&QNF -oo W = * / H 

1 (2n + 1) ± v/1 - 8mV a 2 /£ 2 2a 

(3.132) 

This now has the appropriate limits to reproduce both tanh and sech 2 QNFs. We 
have not found explicit formulae of this type in the extant literature. Note that 
asymptotically 



n 1 ± a/1 - 8mV a 2 /h 2 , nl , N 
fcQNF I± oo(n)-M- + i Y 2q 7 +0(l/n), (3.133) 

in accordance with the general suspicions based on black hole QNMs J3l], ^, |33], |4, 
35| , |36| , |37| ^8| . Finally we consider the quasi- normal energy 



^QNF(n) - l/ +0O + — !/_«, + — , (3.134) 

and compute 

ma 2 (V +OD - V-oo) 2 



E Q Mn) - 2h2 



(2n+l± v/1 - 8mV a 2 /hA 

^4-ocH-V-oo h 2 ^ 2n+1±v / 1 _ 8mV J )a 2 / ^ 2 . (3135) 



2 8ma 2 

Again, we have not found explicit formulae of this type in the extant literature. 

3.10 Related potentials 

(Morse, Poschl— Teller, Manning— Rosen, Hulthen, Teitz, Hua) 

A number of other closely related potentials are also of some interest — though most 
often one is simply revisiting the Eckart potential in disguise. 
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3.10.1 Morse (1929) 

Consider the potential 

V(x) = V (1- exp(-[x - x ]/a)) 2 . (3.136) 

This Morse potential is actually a somewhat odd limit of the (Mobius) 2 potential as 
various parameters go to unity or zero. It is most useful as a model for bound states 
and does not define a scattering problem. 

3.10.2 Poeschl-Teller (1933) 

We should first warn the reader that the actual article by Poschl and Teller [|TS|] is 
somewhat difficult to get hold of. That article starts by discussing the potential 

V(x) = vA , 2 f + 2 f }; xG(0,vra/2), (3.137) 

[sm (x/a) cos z (x/a) ) 



and its hyperbolic analytic continuation (a — > ia 

A B 

sinh (x/a) cosh (x/a 



V(x) = V \ , i2 + — , 2 . — }; x6(0,oo), (3.138) 



relying only on a change of font (typeface) to make the distinction between ordinary 
and hyperbolic trigonometric functions. Finally the article focusses attention on the 
specific case 

V( x ) — 27 — TT = ^o sech 2 (x/a); x G (—00,00). (3.139) 

cosh (x/a) 

Because of this many authors use the phrase "Poschl-Teller potential" to refer to 
the sech 2 potential. While this is historically somewhat inaccurate, insofar as the 
sech 2 potential is already contained as a special case of the Eckart potential, this 
terminology now seems firmly embedded in the literature. (Oddly, Poschl and Teller 
refer to both the Morse 1929 |]IB| and Rosen-Morse 1932 |T7[ articles, but not the 



Eckart 1930 |T§ article.) 

3.10.3 Manning-Rosen (1933) 

Consider the potential 

V(x) = B coth(x/a) - C cosech 2 (x/a); x G (0, 00). (3.140) 

The relevant citation |20 is only an abstract in a report of a conference. To find 



it with online tools such as PRO LA look up Phys. Rev. 44 (1933) 951, and then 
manually scan for abstract #10. The form actually given in the abstract is 

v(x) _ a exp(-2x/&) exp(-s/6) 

V[X} ~ A [1 - exp(-x/6)P + B 1 - exp(-x/b) ' (3 " 141) 
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which one can manipulate into the form above by noting 

2 exp(-ar/6) = l+exp(-g/6) = exp(x/2b) + exp(-x/2b) = _ t w /2 ^ 
1 — exp(— x/6) 1 — exp(— x/b) exp(x/26) — exp(— x/2b) 

(3.142) 

and 

coth 2 ^ = 1 + cosech 2 z. (3.143) 

Note that the Manning-Rosen potential can be obtained from the Eckart potential 
by the formal substitution x — > —x + ma/2 so that 

f = -exp(2x/a) -> +exp(-2x/a). (3.144) 

In particular, Manning-Rosen can be written in the form 

V(x) = A + Bcoth{x/a) + Ccoth 2 (x/a) = A + [B + C coth(x / a)} 2 . (3.145) 

We can get this from the general (Mobius) 2 form of the Eckart potential by appro- 
priately choosing the parameters. Because of the pole at x = the potential is 
best thought of as being defined on (0, oo). It does not define a scattering problem, 
though it may be useful for investigating bound states. 

3.10.4 Hulthen (1942) 

Consider the potential 

/ s exp(— x/a) , , 

V(x) = V Pl ' ; x G 0, oo . 3.146 

1 — exp(— x/a) 



This Hulthen potential |H is actually a special case of the Manning-Rosen potential. 



We can also get this from the general (Mobius) 2 form of the Eckart potential by 
appropriately choosing the parameters. Because of the pole at x = the potential is 
best thought of as being defined on (0, oo). It does not define a scattering problem, 
though it may be useful for investigating bound states. 

3.10.5 Tietz (1963) 



One version of the Tietz potential |22 



is: 



n*) = vJ j^^M V (3 . 147) 

\{smn, cosh, exp} [x/a) J 

We can get this from the general (Mobius) 2 form of the Eckart potential by appro- 
priately choosing the parameters. Depending on the specific choices made it may or 
may not define a scattering problem, though it may be useful for investigating bound 
states. 



-24- 



3.10.6 Hua (1990) 

Hua's potential is |E3 



^) = Vo( 1 l " eXP( 7 2a;/ ?\ V. (3-148) 

\l-qexp{-2x/a)J V ' 



We can get this |2^, |24]] from the general (Mobius) 2 form of the Eckart potential by 



appropriately choosing the parameters. We note 
V(x) = V 

= v n 



exp(x/a) — exp(— x/a) x 
exp(a;/a) — q exp(— x/a) 
sinh(x/a) 



V u 

{sinh, cosh, exp} (x/a) 

(Tietz potential), 

A sinh(x/a) + B cosh(x/a) 



1 + g)sinh(:r/a) + (1 — g)cosh(a;/a) 

(3.149) 

If q > define (1 - q)/(l + q) = tanh^. If q < define (1 + g)/(l - q) = tanh^. 
Then we see 

\{sinh, cosh} (x/a + &)/ 
sinh(x/a — 9) 



i j± sma^x/ a) -t- a cosni^x/ aj \ 
\ {sinh, cosh} (x/a) J 

= (Eckart/ Manning-Rosen as appropriate). 

(3.150) 

Summary: So all of these potentials are either identical to the (Mobius) 2 potential, 
or special cases of the (Mobius) 2 potential. To be historically accurate we should 
really just call this whole collection of potentials the Eckart potential, or appropriate 
special cases of the Eckart potential, as Eckart seems to have been the first author 
to have given the general form. Unfortunately other names are now in such common 
use that historical accuracy is difficult (if not impossible) to recover. 

4. Discussion 

In this paper, we have collected many known analytic results, and described several 
significant new results on analytic QNFs, in a form amenable to comparison with 
the extant literature. In particular we have, in addition to the QNFs themselves, 
focussed on transmission amplitudes, transmission probabilities, and transmission 
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Inter-relationships between various "exactly solvable" potentials 



Name 


Potential V(x) 


Properties 


Morse (1929) 


V (1 - exp(-x/a)) 2 


Special limit of Eckart/(M6bius) 2 


Eckart (1930) 


1 _£ (1-^)2 5 £ = exp(2x/a) 


^ Rosen-Morse ^ (Mobius) 


Rosen-Morse (1932) 


A + 5tanh(x/a) + Csech 2 (x/a) 


«4> Eckart <^ (Mobius) 2 


Morse-Feshbach (1954) 


Vo cosh /i{tanh([x — fia]/a) + tanh/u} 2 


4^ Rosen-Morse 44> Eckart 


Eckart/ Rosen-Morse 


A + Bt&nh(x/a) + Ctanh 2 (x/a) 


4$ (Mobius) 2 function of exp(— 2x/a) 


(Mobius) 2 


V 


' A + B exp(-2x/a)" 
C + Z) exp(-2x/a)_ 


2 


The "best" of these equivalent forms 


Manning-Rosen (1933) 


A + B coth(x/a) — C cosech (x/a) 


Special limit of Eckart/(M6bius) 2 


Hulthen (1942) 


exp(-x/a) 

V ° 1 ( i \ 
1 — exp(— x/a) 


Special case of Manning-Rosen 


Tietz (1963) 


/ sinh([x- x ]/a) \ 2 
\ {sinh, cosh, exp} (x/a) / 


Special limit of Eckart/(M6bius) 2 


Hua (1990) 


y ( l-exp(-2x/a) \ 2 
\ 1 — q exp(— 2x/a) / 


Eckart or Manning-Rosen or Morse 



Table 1: This table shows the inter-connections between many "exactly solvable" poten- 
tials. Many of these potentials are identical to each other, though this may not always be 
obvious at first glance. 



resonances. We did this explicitly for the delta-function potential, double-delta- 
function potential, and asymmetric double-delta-function potential; the step barrier, 
rectangular barrier, and asymmetric rectangular barrier; the tanh potential, sech 2 
potential, and Eckart potential and its variants. In almost all of these cases we 
have been able to take the calculation of the QNFs somewhat further, sometimes 
significantly further, than currently available sources. 

In particular, we have noted that the Eckart/Rosen-Morse/Morse-Feshbach po- 
tentials are actually identical, and that they are generally a (Mobius) 2 function of the 
variable exp(— 2x/a). Indeed many of the "exactly solvable" potentials commonly 
encountered in the literature are actually the same quantity in disguise, typically 
the (Mobius) 2 potential itself or some special case thereof -- and so really should 
just be collectively referred to as variants of the Eckart potential. We should also 
mention that there has recently been some progress in analyzing the (approximate 
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highly damped) QNFs for piecewise Eckart potentials |35], |36], 

What message can we extract concerning the commonly conjectured "(offset) + 
m(gap)" behaviour for highly damped QNFs? From the examples we have seen here 



(and in |35], |36|, |37], p8fl ) is appears that the "(offset) + m(gap)" behaviour depends 
on both a non-zero width for the potential, and a certain amount of smoothness. 
The double-delta potential leads to imaginary gap, imaginary offset, and logarithmic 
sub-leading terms — which is not what one would naively have expected. 

Finally, we reiterate that very few potentials have exact analytically known quasi- 
normal frequencies [QNFs]. Even for so-called "analytically solvable" potentials it is 
not necessarily true that the QNFs can be explicitly located in closed form. Thus 
apart from their intrinsic interest, these exact and approximate results serve as a 
backdrop and a consistency check on ongoing efforts to locate and understand QNFs 
in general physical situations. 
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